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Abstract 

We analyze large sets of energy-release data created by stress-induced brittle fracture in a pure 
sapphire crystal at close to zero temperature where stochastic fluctuations are minimal. The 
waiting-time distribution follows that observed for fracture in rock and for earthquakes. Despite 
strong time correlations of the events and the presence of large-event precursors, simple prediction 
algorithms only succeed in a very weak probabilistic sense. We also discuss prospects for further 
cryogenic experiments reaching close to single-bond sensitivity and able to investigate the existence 
of a transition-stress regime. 

PACS numbers: 62.20.Mk,91.30.Px,07.20.Mc 
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We have recently described a serendipitous and novel measurement of brittle fracture 
using cryogenic calorimetryi. In a 260 g pure sapphire crystal cooled to 20 mK, cracks formed 
under pressure from sapphire bearings in what amounts to a sharp indentation experiment 
(Fig.[Tl[2]). The small contact surfaces generated stress fields vanishing quickly with distance, 
and ensured stability of the fractures^. The calorimetric measurement provided a direct 
measurement of the energy of the phonons from fracture events, and great sensitivity, of 
the order of a few femto- Joules. The rich and complete event catalogues, of many thousand 
femto-fractures each, contain the arrival time and energy of each event, and show several 
statistical similarities to earthquakes, despite the many orders of magnitude difference in the 
energy ranges^. The similarities include: (i) the probability distribution of fracture-energy 
release is a power law with an exponent close to that of the differential Gutenberg-Richter 
relation expressed for seismic moment (which is proportional to earthquake energy)^, [ii) 
fracture events are long-range correlated in time with a power-law waiting-time distribution 
for short times, {in) the fracture time series has the characteristics of fractal Gaussian 
intermittant noise, and {iv) there is an elevated event rate right after large events and a 
power-law event rate decay. More generally, the absence of trends in the data indicate that 
this represents a new example of steady-state slow brittle fracture, in an ordered system. Up 
to now, such fracture has been linked to the disorder inherent in self-organized-criticality^i^i^. 

In the following, we show that the waiting-time distribution follows a general power-law 
exponential form observed in earthquakes and rock fracture, with the same power. We 
demonstrate further correlations in the data and attempt to use them as predictors of the 
large, catastrophic fractures that should eventually occur. Lastly, we discuss a dedicated 
cryogenic experiment to study these phenomena, down to energies close to those of single 
bonds in the crystal. Such an experiment could also probe putative stress-dependent varia- 
tions of the fracture rate, and investigate the existence of a transition-stress regime where 
the average energy release would vanish. 

I. WAITING TIMES AND CLUSTERING 

We have showni"^ that the distribution of the waiting time between consecutive events 
above threshold, w, follows a power law at short times with an exponential fall-off at large 
waits: dN/dw oc w'"' exp — w/wq. The form of this waiting-time distribution is identical 
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to that observed for earthquakes and rock fracture, further extending its vahdity^. Similar 
forms may be derived from the Gutenberg-Richter and Omori laws^. The average wait 
must be proportional to the scale term wq; integration by parts yields: w = {1 — a)wo. 
On the other hand, the distribution of events as a function of energy follows a power law: 
dN/dE oc E~l^ , with (3 ~ 1.9^. By integration, the number of events above a given energy 
therefore also follows a power law: A^(> E) oc E~^^^. This is inversely proportional to 
the average waiting time for events above a threshold: w{> E) = {1 — a)wo{> E) oc 
Fig. [3] shows that, as threshold energy increases, the distribution of waiting times retains 
the form exp {—w/wo{> E)), where the power a has little dependence on the threshold 
energy and Wq scales like E^~^. A fit of the wait power yields a = 0.33 ± 0.01. With the 
notations from Ref.^, we find 5 = = 1.49 ± 0.02 and 7 = 1 - a = 0.67 ± 0.01. These 
values are strikingly close to those obtained for earthquakes^ and are also compatible with 
those for acoustic emission of fractures in rocks^. 

If the energy distribution can be extrapolated to values large enough to cause a catas- 
trophic destruction of the detector itself, then such a catastrophic event would arrive in 
a long, perhaps, but finite, time. If Ecat is the energy released as the crystal breaks, 
then, in this setup, it would be expected after = Wo{> E) [Ecat/ E]^"^ ~ (0.003 h) x 
[£'cat/(10 keV)]°'^. This scales slightly slower than a linear relation. To obtain an order of 
magnitude of the timescales involved, we assume the fracture surface energy of sapphire is 
an upper limit on the energy that would be released by a crack. For instance, taking a value 
of 7.3 J/m^ ^ 4.5 x 10^^ keV/cm^ for the surface energy of the {1012} plane^i, and defining 
a catastrophic crack size as 1 cm^ in the 4x4x4 cm^ cubic crystal, this translates to an 
upper bound of E^at < 4.5 x 10^^ keV. The weak upper limit on the expected wait for such 
a catastrophic event is therefore several millenia. 

The distribution of time intervals between all, rather than consecutive, events above 
various energies is shown in Fig. HI As energies increase, the distribution becomes more and 
more peaked at low time intervals. This is further indication that the large events cluster-. 
As a control, the same analysis is applied to a random shuffle of the arrival times in the 
data. In the shuffle, there is no energy- dependent effect. 



3 



II. WEAK PREDICTABILITY 



Some examples of event-energy time series from a 50 h run containing ^ 30000 events 
above a threshold of 12 keV are shown in Fig. O Various cases are visible, including a large 
event with a precursor {t ~ 16.154 h), a relatively isolated large event (t ~ 22.915 h), and a 
lull before a large event followed by aftershocks (t ~ 22.93 h). In Figures[6]A and A', we plot 
the average value of the waiting time before each of the ~ 14000 small events (12-30 keV) 
and ~ 1100 large events (300-1000 keV). The waiting times for small (respectively large) 
events are here defined as the wait between a small (resp. large) event and the preceeding 
event regardless of its size. On average, there is less wait before large events (0.0013 h) 
than before small events (0.002 h). To check if this a statistical fluctuation, we generate 100 
shuffles of the data set, by randomly permuting the arrival times of the events, then do the 
same analysis as on the original data. The distribution of the average values of the shuffles 
does not cover the spread of the real values, confirming that in the actual data, the wait is 
shorter before large events than before small ones. This appears to be another manifestation 
of the increase in rate which is particularly evident around large events^. 

These and the numerous other correlations present in the data provide motivation to 
attempt prediction of large events, a challenge of relevance for other phenomena, ranging 
from avalanches in snow^^ to earthquakes"^. Fig. [6] B compares the distribution of waiting 
times before small and large events. The significant correlations present on average are much 
harder to exploit on an event-per-event basis, as the distribution for large events does not 
differ greatly from that for small events. 

We also attempt to predict the arrival of large events using the distribution of events in 
a given time window (Fig. [6]C). Window duration is 0.002 h, corresponding to the average 
waiting time in the run. For comparison, we generate 5000 random intervals. The difference 
between the distribution of counts in the random intervals and in the intervals preceeding 
small or large events is slim, while the difference between intervals preceeding small and 
large events is slighter yet. With these simple methods, predictability of individual large 
events is therefore poor. 

While the weak predictability we have described here could perhaps be enhanced by more 
sophisticated algorithms, it might also be either a general conclusion for brittle fracture, or 
indicate that, in our particular setup, the combination of several crack systems propagating 
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independently masks any individual patterns and predictability. 

III. PROSPECTS FOR FURTHER STUDY 

Further study would benefit from a dedicated cryogenic detector with only a single bearing 
creating the cracks. It could allow investigation of fractures down to low energies. Addi- 
tionally, if fracture rate is found to depend on applied stress, it might allow investigation of 
low rates close to a possible transition stress at which fractures just appear. If such a regime 
exists, the waiting-time distribution would be a pure power-law, and the whole system could 
be in a critical transition at which the average energy release rate vanishes. 

One option would be to carry on in a low-background, deep-underground setup, such 
as that of the CRESST II experimenti^. An existing detector holder could be modified to 
include a single bearing pressing against one end of the cylindrical crystal of 40 mm height 
and 40 mm diameter. In itself, this will require some ingenuity as the crystal itself must not 
move because of the bearing but cannot, for thermal reasons, be held firmly by large contact 
areas. One of the currently standard CaW04 crystals could be used, or a new AI2O3 one 
could be manufactured. In either case, it would be interesting to retain the light detector 
of the CRESST II setup to see if crack formation is accompanied by light emission, since 
fracto-emission of photons and electrons has been reported in other crystals^"^. Adjusting the 
tightness of the spring pressing the bearing to probe an effect on crack rate, and to reach 
transition stress if it exists, would have to be done between cryogenic cycles and would 
require some trial and error. As in the original work, energy calibration would be obtained 
by an external, removable, ^''Co source, providing 122 keV photons. Lower-energy calibration 
would be obtained via heater pulses. It should be possible to lower the phonon threshold 
down to around 1 keV compared to the ~ 10 keV in this work, though this gain would 
probably not be significant from the standpoint of brittle fracture. Another underground 
setup that could be of interest is that of the EDELWEISS experiment which uses germanium 
ionization-phonon detectors^-. In this case, there would be a simultaneous measurement of 
the phonons created by the cracks as well as of whatever ionization the cracks create^. 

Though the original work was carried out in a special low-background environment un- 
derground, it would be simpler if in the future it could be done in a standard cryostat on the 
surface. For this, the rate of crack events must be much larger than the rate of other events; 
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this means that though a low-threshold experiment may be feasible, a transition-stress one 
will be difficult. The rate of crack events may depend on the force pressing the bearing 
against the crystal but does not depend on the mass of the crystal itself. We assume that 
in a dedicated experiment, there would be only a single bearing rather than a dozen as in 
the original work, and therefore divide the original rate by 12. The competing backgrounds, 
mainly cosmic-ray-induced particles and radioactivity of the detector and its surroundings, 
both increase with detector mass. As illustration, we extract a rough estimate of the sum 
of these backgrounds from previously published data for a partially shielded 1 kg sapphire 
detector at the Earth's surface^. Between the threshold of 50 keV and ~ 1 MeV, the back- 
ground follows the product of an inverse power law, with an exponent of 0.6, and a decaying 
exponential, with a typical energy of ~ 600 keV, high above the range relevant here. In 
Fig. [3, we compare the scaled rate of cracks to the backgrounds in an 0.1 g sapphire detec- 
tor for two different scaling laws of the background as a function of detector mass m: scaling 
proportional to mass (i.e. volume, oc m), and scaling proportional to surface area (oc m?^^). 
In both cases, we assume that the background power-law holds below 50 keV. A significant 
background proportional to the surface area comes from cosmic muons of which there are 
^ 1 cm~^ min~^ depositing ^ 500 keV per mm of AI2O3 passed through^. Working at 
ground-level will require some combination of small crystals, low threshold, and, if possible, 
increased crack rate, though this last point is incompatible with a reduction in stress to the 
level at which fractures may nearly vanish. Energy calibration of a small crystal requires a 
low-energy radioactive source such as ^^Fe; the source must be placed within the cryostat 
for calibration, though it may be possible to remove it or block it during data-taking for 
fractures. Such a small crystal could have a threshold of less than a few hundred eV, and 
as energy decreases, the crack rate (oc E~^'^) rises faster than the extrapolated background 
(oc 

In addition, an energy threshold below a few eV could make the device sensitive to the 
rupture of single sapphire bonds. A binding energy of 7.34 eV per atom has been reported for 
A\.20^, though it is not apparent to us what phonon energy accompanies rupture of a bond. 
A 10 eV threshold is achievable given current, ~ cm^, cryogenic detector development. For 
instance, the CRESST experiment has developed thin silicon calorimeters, of surface area 
several cm^ and thickness about 0.5 mm, with thresholds better than 40 eV— These 
detectors are optimized for light detection rather than for a low threshold per se. A smaller. 
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parallelipiped-shaped, AI2O3 device, measuring 5x5x1 mm^ (m=0.1 g), with an optimized 
thermometer, should be able to reach lower thresholds, while remaining large enough for a 
pressure-bearing 1 mm diameter sapphire or diamond ball. The device could be mounted 
^ 1 mm from a CRESST-type light detector of similar size, to see any light produced by the 
fractureslA. We note that to obtain an absolute energy calibration at 6 keV from ^^Fe and to 
have a threshold of 10 eV will require a dynamic range of about three orders of magnitude 
which may be difficult to obtain with a transition-edge sensor. Another challenge will come 
from the rate of cracks. Detectors of size 3 x 3 x 0.5 mm^, optimized for speed rather than 
threshold, reach rise times of ~ 1 /is^^ (smaller devices in which the transition-edge sensor 
itself is the absorber can be an order of magnitude faster—) . Though this may be compatible 
with the average rate of cracks, pileup will be inevitable, even for an arbitrarily fast detector, 
since the distribution of waiting times contains an inverse-power-law term (Fig. [7]). 

The calorimetric technique is readily applicable to many other dielectric materials, such 
as CaW04, Ge and Si already mentioned. As we have stated, we are not however aware of 
a clear relationship between the elastic energy used to break bonds and the elastic energy 
left over in phonons which we measure. Nonetheless, the partition of energy is simpler 
than in the case of acoustic emission, where only a fraction of phonons are measured. The 
calorimetric technique could therefore provide new insight into the mechanics of fracture. 

IV. CONCLUSION 

The distribution of waiting times between brittle fracture events observed in a cryogenic 
detector contains a power-law term which is independent of energy threshold and an expo- 
nential scale that depends on it. This form matches that previously observed by acoustic 
emission in rock and that observed for earthquakes. If the energy distribution holds for 
events large enough to shatter the detector, then such an event is expected in a very long, 
but finite amount of time. Though we have shown additional correlations in the data, pre- 
dicting such large, catastrophic, fractures is not straightforward. To see if this is due to 
the multiple sources of cracks in this data, we propose a dedicated experiment with a single 
pressure point. With typical fracture rates observed heretofore, such an experiment is fea- 
sible at ground level with a smaller cryogenic detector if extrapolation of background holds 
to low energies and masses. However, searching for vanishing fracture rates requires at least 



7 



a shallow underground site. In either case, the lower threshold associated with a smaller, 
optimized, detector would enable it to probe brittle fracture down close to the energy of 
single bonds in the crystal. The calorimetric technique could provide additional insight into 
the partition of elastic energy into permanent dislocations and phonons. 
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FIG. 1: Confocal microscopy picture of sapphire crystal fractured by sapphire bearing. Diame- 
ter of affected area is pa 2 mm. Shde marks are visible, as are irregular fractures of radial and 
circumferential type. 
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FIG. 2: Fractured area of the crystal scanned by confocal microscopy to a depth of about 150 /im. 
Features visible at the sapphire surface are openings of two-dimensional cracks that extend fairly 
deep into the crystal. Their depth could not be measured accurately with the present method. This 
three-dimensional image was made by volume rendering of the cracks only, i.e. the solid material 
is not shown here. 
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FIG. 3: Left: distributions of waiting times w between consecutive events above various energy 
thresholds. Error bars are the square root of each bin content. Data are well fitted by the product 
of an inverse power law and an exponential decay (oc w^°^ exp —w/wq). Top right: the wait 
distribution scale term, wo{> E), obtained from fits of the wait distribution (error bars are from 
fit) is compatible with fit by oc E^~^ (dashed red line), where (3 = 1.9 is obtained from the fit of the 
energy distributior>i. Bottom right: wait exponent a obtained from fits, as a function of threshold 
energy (error bars are fit errors). Data are fitted by constant function yielding a = 0.33 it 0.01. 
These figures demonstrate that at least up to the highest energies, the wait scale term does indeed 
scale like wo{> E) oc El^~^ , whereas the wait power term does not depend strongly on the energy 
threshold. 12 



Data 



■ -I- - I 



Counts per Unit Time (AU) 

0.16T-- 
0.13 
0.11 
0.09 

0.07 
0.06 

0.05 
0.04 



E>12 keV 
E>30 keV 
E>100 keV 
E>300 keV 




H 1 — \ — I \ \ \ \ 1 \ \ — I — I 1 1 1 1 \ H 

0.001 0.01 

Time Interval Between All Events Above Threshold (h) 



0.1 



Shuffled Data 



Counts 

0.16 

0.13 
0.11 
0.09 

0.07 
0.06 

0.05 
0.04 



per Unit Time (AU) 



■ -I- -i 





H \ \ — I — |- 



H 1 1 1 \ h 



0.001 0.01 

Time Interval Between All Events Above Threshold (h) 



0.1 



FIG. 4: Distribution of times between all (as opposed to consecutive) events in various energy 
ranges (note that binning varies between ranges). Top: in the data, the distribution becomes 
peaked at low times as the event size increases. This is an indicator of clustering for large events. 
Bottom: in a random shuffle of the data, the distributions differ less for all event sizes. 
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FIG. 5: Examples of event energies as a function of time. Some large events {E > 300 keV) appear 
after relatively quiet periods (for instance just before 22.92 h) whereas others have some precursors 
(for instance just before 16.155 h) and others display aftershocks (around 22.93 h) 
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A) Before 12-30 keV Events 
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FIG. 6: Fig. A shows the average wait between a small (12-30 keV) event and the event of any 
size preceeding it (solid blue bar), and the distribution of this average wait before small events in 
100 shuffles of data (dashed red histogram). Fig. A' shows the same but for large events (300- 
1000 keV). Fig. A and A' demonstrate that there is significantly less wait on average before a large 
event than before a small one. However, Fig. B illustrates that though the averages differ, the 
distribution of waiting times before small and large events are quite similar and will not provide 
strong discrimination between individual small and large events. Fig. C shows the distribution 
of counts in random intervals (solid blue line), in intervals before small events (dashed red) and 
in intervals before large events (dotted green). The distributions do not show a strong distinction 
between small and large events, and indeed littlleBdifference between random intervals and intervals 
preceeding events. 
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FIG. 7: Top: crack rate compared to backgrounds in hypothetical 0.1 g detector in the case of a 
background scahng hke detector mass and in the case of one scahng hke detector area. Rate of 
cracks is scaled to a single pressure point, though in a new experiment this rate could perhaps 
be fine-tuned between cryogenic runs by altering the force applied on the bearing. Backgrounds 
are natural radioactivity and cosmic-induced at ground level, and are extrapolated from Fig. 2 of 
Ref.^^ to low energies and to low detector masses. Being able to observe the cracks at ground level 
will require a small crystal and perhaps an increased rate of cracks. Bottom left: the crack rate 
integrated above even a low threshold remains, on average, compatible with typical thermal time 
constants of cryogenic detectors and typical acquisition rates. Bottom right: the distribution of 
waiting times, here in the oc w"" transition-stiidss limit with a = 0.33, shows that the abundance 
of short waiting times makes pile-up inevitable. 



